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Analysis  of  a  signal  starting  time 
estimator  based  on  the  Page  test 
statistic 

D.  A.  Abraham 


Executive  Summary:  In  the  detection,  localization,  and  classification  of 
targets  in  anti-submarine  warfare  using  active  sonar,  the  presence  of  a  target 
echo  must  first  be  detected  in  the  received  array  data.  From  the  information 
obtained  by  the  detection,  the  target  is  marginally  localized  in  range  by  the 
time  of  detection  and  in  bearing  by  implementing  detectors  on  multiple  beams 
that  selectively  accept  energy  from  specific  directions.  The  target  echo  must 
then  be  segmented  in  the  received  array  data  prior  to  refining  the  initial  bearing 
and  range  estimates,  depth  estimation,  or  target  classification.  This  requires 
estimation  of  the  start  and  stop  time  of  the  target  echo,  a  particularly  difficult 
problem  when  the  effect  on  the  transmitted  signal  due  to  propagation  through 
the  ocean  medium  and  reflection  off  the  target  are  unknown  prior  to  localization 
and  classification. 

In  this  memorandum  a  starting  (or  stopping)  time  estimator  based  on  the  Page 
test  statistic  is  proposed  and  analyzed.  The  Page  test  is  a  sequential  test  used 
to  detect  a  change;  in  this  case  the  transition  from  signal- absence  to  signal- 
presence.  The  effects  of  signal  power  and  the  Page  test  detection  threshold  are 
investigated  for  two  signal  models.  The  results  of  this  analysis  may  be  used 
to  ensure,  within  some  probability,  that  the  target  echo  segmentation  includes 
the  complete  echo. 
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Analysis  of  a  signal  starting  time 
estimator  based  on  the  Page  test 
statistic 

D.  A.  Abraham 


Abstract:  The  time  of  the  most  recent  reset  to  zero  of  the  Page  test  statis¬ 

tic  is  proposed  as  an  estimator  of  the  starting  time  of  a  signal.  The  probability 
mass  function  of  the  estimator  is  determined  analytically  subject  to  a  quanti¬ 
zation  of  the  Page  test  statistic  update.  Closed  form  results  for  the  first  three 
uncorrected  moments  of  the  estimator  are  presented.  The  analytical  results  are 
verified  by  comparison  with  simulation  results  and  the  fineness  of  the  quanti¬ 
zation  required  for  accurate  representation  is  investigated  by  evaluation  of  the 
Kolmogorov- Smirnov  statistic.  The  bias,  standard  deviation,  and  skewness  of 
the  estimator  as  a  function  of  the  signal  strength  and  detector  threshold  are 
evaluated  for  Gaussian  shift-in-mean  and  noncentral  chi-squared  signal  types. 

Keywords:  Page  test  o  submarine  detection,  localization,  classification  o  sig¬ 
nal  starting  time  estimator 
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1 

Introduction 


Many  signal  processing  applications  require  the  detection  of  a  change  occurring  in 
a  sequence  of  observed  data  [1].  One  type  of  application  may  be  represented  by  a 
model  where  the  time  of  occurrence  of  the  signal  is  unknown  and  once  the  signal 
starts,  it  does  not  stop  until  appropriate  action  has  been  taken.  Examples  of  this 
type  include  the  detection  of  mechanical  failures  in  machinery,  the  detection  of 
changes  in  heart  sounds,  and  the  rapid  detection  of  approaching  ships  or  torpedoes 
in  sonar  applications.  Another  type  of  change  detection  may  be  represented  by  a 
signal  model  where  the  duration  of  the  change  is  finite  and  both  the  duration  and 
time  of  occurrence  are  unknown.  Examples  of  this  type  include  the  segmentation 
of  words  in  speech  processing  and  the  detection  of  active  sonar  echoes  in  shallow 
water.  In  most  of  these  examples,  once  the  signal  is  detected  some  action  is  taken 
that  requires  an  estimate  of  when  the  signal  started  (and  possibly  when  it  has 
stopped).  For  instance,  the  components  produced  by  a  machine  may  need  to  be 
examined  for  adherence  to  design  specification  once  a  failure  has  occurred,  a  word 
may  need  to  be  classified  once  it  has  been  segmented  in  speech,  or  an  active  sonar 
signal  echo  may  be  localized  in  depth  or  classified  once  it  has  been  detected  and 
segmented. 

To  mathematically  introduce  the  change  detection  problem,  suppose  that  the  data 
Xi,  X2,  ...  are  observed.  For  time  indices  i  <  ns  assume  the  data  is  distributed 
according  to  a  null  or  before-change  rule  ( Hq )  and  for  i  >  ns  assume  the  data 
is  distributed  according  to  an  alternative  or  after-change  rule  (#1).  It  is  desired 
to  determine  when  the  change  has  occurred  with  as  little  data  as  possible  and  to 
subsequently  provide  an  estimate  of  the  time  of  change,  ns. 

The  Page  test  [2]  has  long  been  used  to  determine  when  the  distribution  of  a  sequence 
of  random  variables  has  changed.  The  recursive  form  of  the  test  declares  that  a 
change  has  occurred  when  the  statistic 

Wk  =  m^{OyWk^+g(Xk)}  (1) 

crosses  a  threshold  where  Wq  =  0  and  g  (•)  is  the  Page  test  update  or  nonlinearity. 
The  Page  test  may  also  be  viewed  as  consecutive  Wald  sequential  probability  ratio 
tests  (SPRTs)  terminating  at  Hq  followed  by  one  terminating  at  H\.  From  this 
interpretation  of  the  Page  test,  a  natural  estimator  of  the  starting  time  of  the  signal 
is  the  time  of  the  most  recent  SPRT  termination  at  Hq;  that  is,  the  time  of  the 
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most  recent  reset  of  W k  to  zero.  Mathematically  the  starting  time  estimator  may 
be  represented  by 

M  =  max  {k  <  N  :  Wk  =  0}  ,  (2) 

where  N  is  the  stopping  time  of  the  Page  test, 

N  =  inf  {&  >  0  :  Wk  >  h}  .  (3) 

A  direct  application  of  the  maximum  likelihood  estimation  principle  to  this  problem 
is  not  possible  unless  the  change  has  already  been  detected  resulting  in  a  fixed  sample 
size  estimation  problem.  Under  this  condition,  the  proposed  estimator  is  shown  to 
be  related  to  the  maximum  likelihood  estimator  (MLE).  First,  define  the  cumulative 
summation  of  the  Page  test  update  as 

k 

Sk  =  J2d(Xi).  (4) 

i=  1 

Now  the  Page  test  may  be  described  non-recursively  as 

wk  =  Sk  —  min  Si 

i<k 

k 

=  E  9(Xj).  (5) 

-  j=i+l 

From  this  description  of  the  Page  test,  the  proposed  starting  time  estimator  is 

N 

M  =  arg  max  ^  g(Xj) .  (6) 

-  j=i+ 1 

If  the  Page  test  update  g  (x)  is  the  log-likelihood  ratio,  then  Eq.  (6)  has  the  same 
form  as  the  MLE  of  the  signal  starting  time  except  that  here  the  sample  size  N  is 
a  random  variable  determined  by  the  termination  of  the  Page  test. 

For  the  detection  and  segmentation  of  finite  duration  signals  one  may  use  consecutive 
Page  tests  with  alternating  hypotheses  [3].  In  this  case,  the  stopping  time  of  the 
signal  is  similarly  estimated  using  the  time  of  the  most  recent  SPRT  termination  at 

This  paper  develops  a  method  for  determining  the  probability  mass  function  (PMF) 
of  the  above  estimator  by  first  quantizing  the  Page  test  update  g  ( x )  as  shown  in 
Sect.  2.  Equations  are  presented  for  the  first  three  moments  from  which  the  variance 
and  skewness  may  be  determined.  In  Sect.  3  the  analytical  results  are  verified  by 
comparing  the  PMF  and  cumulative  distribution  function  (CDF)  with  simulation 
results.  Additionally,  the  Kolmogorov-Smirnov  statistic  is  determined  as  a  function 
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of  the  fineness  of  the  quantization  of  the  Page  test  update.  Finally,  in  Sect.  4  the 
performance  of  the  estimator  is  determined  as  a  function  of  the  threshold  used  in 
the  Page  test  and  the  signal  strength  for  Gaussian  shift-in-mean  and  noncentral 
chi-squared  signal  types. 
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2 

Estimator  analysis 


2.1  Quantization 

Following  the  method  of  Brook  [4]  and  as  further  described  in  reference  [5],  the 
Page  test  update  g  ( x )  is  quantized  into  levels  equally  spaced  at  intervals  of  width  A. 
Arbitrarily,  let  the  levels  be  k  =  iA  for  i  -  0,  ±1,  ±2,  . . ..  Let  the  quantization  rule 
be  2  ( 9 )  =  k  if  g  £  [A  (i  -  0.5) ,  A  (i  +  0.5)),  where  Q  (•)  represents  the  quantization 
operation.  Let  the  integer  7  be  such  that  77_i  <  h  <  where  h  is  the  threshold  for 
the  Page  test.  Thus,  while  the  Page  test  is  running  (before  a  threshold  crossing), 
the  quantized  statistic  must  lie  in  one  of  the  states  {/ij/Jo- 

Suppose  that  under  Hj  for  j  =  0  or  1  the  Page  test  update  has  probability  distri¬ 
bution  function  (PDF)  fj  ( g )  and  CDF  Fj  (g).  Then  the  probabilities  of  observing 
the  levels  under  Hq  are 

Pi  =  Pr0  (A  (i- 0.5)  <g<  A  (*  +  0.5)} 

=  F0  (A  (i  +  0.5))  -  F0  (A  (i  -  0.5))  (7) 

and  under  Hi  are 

Qi  -  Prj  (A(*-0.5)  <g  <  A(i  +  0.5)} 

=  Fi  (A  (i  +  0.5))  —  F\  (A  (i  —  0.5)) .  (8) 


2.2  Probability  mass  function 

There  are  two  possible  scenarios  that  must  be  considered  separately  in  determining 
the  PMF  of  the  starting  time  estimator  M :  (i)  when  the  Page  test  statistic  does 
not  reset  to  zero  before  crossing  the  threshold  once  the  signal  starts,  and  (ii)  when 
it  does  reset  to  zero  before  crossing  the  threshold.  Without  a  loss  of  generality  the 
actual  signal  start  time  can  be  subtracted  from  all  time  indices  and  the  starting 
time  estimator  M  so  only  the  estimator  error  is  considered.  Case  (i)  then  represents 
negative  values  of  M  and  case  (ii)  positive  values  of  M.  Note  that  the  data  now 
follows  Hq  for  times  k  <  0  and  follows  77 1  for  times  A:  >  0. 

Let  the  variable  u  represent  a  vector  of  length  7  composed  of  the  probability  of 
being  in  states  lo,  li,  ■  ■  ■ ,  Z7-i-  Let  the  variable  v  represent  a  vector  of  length  7  —  1 
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composed  of  the  probability  of  being  in  states  l\)  I2,  . . . ,  /7_i*  These  variables  will 
be  appended  with  subscripts  to  distinguish  various  conditions  or  time  dependence. 
The  u  variables  will  be  referred  to  as  state  probability  vectors  and  the  v  variables 
as  reduced  state  probability  vectors.  Also  note  that  they  may  be  related  by  the 
matrix-vector  product 

V  =  [0  Iy_l]u 

=  PU,  (9) 

where  P  is  a  7  -  1  by  7  matrix  composed  of  a  length  7  —  1  vector  of  zeros  (0)  and 
a  dimension  7  —  1  identity  matrix  (I7_i). 


Let  S  (uo)  represent  the  event  that  the  Page  test  statistic  proceeds  directly  to  the 
threshold  without  entering  the  state  Zo  when  the  state  probability  vector  is  initially 
uo  and  the  data  follows  condition  H\  (i.e.,  the  change  has  occurred).  In  this  situ¬ 
ation,  state  lo  may  only  be  observed  at  the  initial  time  sample,  say  k  =  0.  Thus, 
for  all  times  k  >  1,  the  reduced  state  probability  vector  may  be  related  to  the  most 
previous  one  by 

vfc  =  DiVfc_i,  (10) 


where  is  the  reduced  state  probability  vector  for  time  sample  k  and  Di  is  a  7  -  1 
by  7  —  1  matrix  of  the  probability  of  a  transition  from  one  state  to  another  under 

Hu 


Qo 

9-1  • 

'  •  9-7+2 

Qi 

qo  ■ 

9-7+3 

?7-2 

?7-3 

9o 

(ii) 


The  length  7  state  probability  vector  for  time  k  =  0,  uo,  may  be  partitioned  into 


uo  = 


(12) 


where  ua  =  e^uo  is  the  probability  of  being  in  state  Zo,  eo  =  [1  0  •  •  *  0]T  is  a 
dimension  7  or  7  —  1  unit  vector  pointing  toward  the  first  coordinate,  and  vo  =  Puo. 
Then,  the  reduced  state  probability  vector  for  time  k  =  1  is 


vi  =  uad  1  4-  Div0 
=  [di  Di]  u0 

=  Diu0,  (13) 

where  Di  =  [di  Di],  di  is  a  vector  composed  of  the  probability  of  a  transition  from 
state  lo  to  the  states  l\,  . . . ,  Z7_i  under  Hu 

di  =  [qi  q2  ••  Q'r-if  ■  (14) 
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Now  the  probability  of  a  straight  climb  to  the  threshold,  without  entering  state  Zo> 
is  the  probability  of  reaching  or  exceeding  the  threshold  on  the  next  update  from 
each  state  times  the  probability  of  being  in  each  state  for  each  future  time  sample, 


Pr{5(u0)} 


where  ca  =  <?;, 


and 


w  = 


oo 

=  Caua  +  CTV0  +  ^2  C'^Vk 
fc=  1 

OO 

=  caua  +  cTv0  4-  cT  ^  Dj-1vi 
fc= l 

=  caua  +  cTv0  +  ct(I7_i  -  Di)_1  Diuo 

T 

=  w  u0, 


SS7-I  Qi 
EiS  7-2 


E  2 


£1  Qt  J 


+  Df(l7-i-Di)"Tc- 


(15) 


(16) 


(17) 


Note  that  the  development  of  Eq.  (15)  requires  that  the  eigenvalue  with  the  max¬ 
imum  absolute  value  of  Di  be  less  than  one.  This  property  will  be  established  in 
Sect.  2.3. 


Now  consider  case  (i)  where,  once  the  signal  starts,  the  Page  test  statistic  does  not 
reset  to  zero  before  crossing  the  threshold.  Let  Am  for  m  <  0  represent  the  event 
that  Wm  =  0  and  that  Wm+ 1,  Wm+ 2,  . . . ,  Wo  are  not  equal  to  zero.  Then  for  m  <  0 
the  PMF  of  M  is 


Pm  M  =  Pr  {S  (no)  |Ai }  ,  (18) 

where  u0  —  [(')  v(5'J  is  the  state  probability  vector  at  time  k  =  0.  The  reduced  state 
probability  vector  at  time  k  =  0  (vo)  is  now  determined  subject  to  the  occurrence 
of  event  *4m.  The  transition  of  the  states  from  time  k  —  ltofcform  +  l<fc<0 
may  be  represented  by 

vfc  =  D0vfc_j,  (19) 


where  Do  is  a  7  —  1  by  7  —  1  matrix  of  the  probability  of  a  transition  from  one  state 
to  another  under  Ho , 


Do 


Po 

V- 1  * 

•  •  P- 7+2 

Pi 

po  • 

••  P- 7+3 

P'1-2 

Vl~  3  • 

Po 

(20) 
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Since  event  Am  requires  that  Wm  =  0  and  Wm+\  ^  0,  the  transition  of  the  states 
from  time  to  to  time  m  + 1  is  restricted  to  moving  from  state  Iq  to  the  states  l\,  . . . , 
Z7_ i  under  Ho.  Thus,  the  reduced  state  probability  vector  for  time  k  =  m  +  1  is 

vm+1  =  (e%uss)  d0,  (21) 


where 


d0  =  [pi  p2  ■■■  p7- i]T  ,  (22) 

eguss  is  the  steady  state  probability  of  observing  state  lo,  and  uss  is  the  steady 
state  probability  of  observing  the  states.  As  derived  in  reference  [5],  uss  is  the 
eigenvector  (normalized  to  sum  to  one)  associated  with  the  maximum  eigenvalue  of 
the  probability  transition  matrix  including  regulation  at  zero  under  Ho, 


St=— oo  Pi 

^2i=-ooPi 

Zji=-oo  Pi 

Pi 

Po 

P- 7+2 

C0  = 

P2 

Pi 

P— 7+3 

(23) 

.  P',-1 

P',-2 

Po 

Now,  the  reduced  state  probability  vector  for  time  k  —  0  subject  to  the  occurrence 
of  event  Am  may  be  determined  by  combining  Eqs.  (19)  and  (21), 

v0  =  (e^uss)  Dq  m_1d0.  (24) 


Now  consider  case  (ii)  where,  once  the  signal  starts,  the  Page  test  statistic  resets 
to  zero  before  crossing  the  threshold.  Let  Bm  for  to  >  0  represent  the  event  that 
Wm  =  0  and  that  at  time  k  =  0  the  Page  test  statistic  is  in  a  steady  state  condition 
(i.e.,  the  state  probability  vector  is  uss).  Then  for  m  >  0  the  PMF  of  M  is 

Pm  M  =  Pr  {<$  (e0)}  Pr  {J Bm }  .  (25) 

Let  ii*.  be  the  state  probability  vectors  for  k  =  0,  1,  ....  m.  As  described  in 
reference  [5],  they  are  related  by  the  product 


u*;  =  CjUfc_i,  (26) 

where  C\  is  the  probability  transition  matrix  including  regulation  at  zero  under  Hi, 


52»=-oo  Qi 

1  — — - —  Qi 

Qi 

9o 

9-7+2 

Ci  = 

92 

9i 

9-7+3 

97-1 

97-2 

9o 
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The  probability  of  observing  state  lo  at  time  m  when  the  Page  test  is  in  its  steady 
state  at  time  k  =  0  is  then 

Pr{Bm}  =  ejum 

=  eZC?uss.  (28) 


Combining  Eqs.  (15),  (18),  (24),  (25)  and  (28)  results  in  the  following  definition  for 
the  probability  mass  function  of  the  error  of  the  stopping  time  estimator  of  Eq.  (2), 


r  /  (eo  u$s)  wrPTD0 

PM  [m]  —  \  f  T  \  Tf^rn 

[  (eo  wj  Bq  C?us 


— ra— 1 


m  <  0 
m  >  0 


(29) 


2.3  First  three  moments 

Define  the  uncorrected  moment  of  the  error  of  the  starting  time  estimator  as 
M*  =  E[M*\ 

=  (-1)*  (e^uss)  wTPT  (  J2  ra'D™-1)  d0 

Vm— 1  / 

+  (eo w)  eo  ^ct]  uss-  (30) 

\m=0  / 

Evaluation  of  Eq.  (30)  for  i  =  0,  1,  25  and  3  is  straightforward  using  closed  forms 
for  the  infinite  summations  as  found  in  reference  [6].  Note  that  this  requires  that 
the  eigenvalue  with  the  maximum  absolute  value  of  the  matrices  Do  and  Ci  be 
less  than  one.  With  the  use  of  Perron’s  theorem  [7],  which  dictates  that  all  the 
eigenvalues  lie  on  a  disk  of  radius  equal  to  the  maximum  absolute  column  sum,  it 
suffices  to  note  that  all  the  entries  in  the  matrices  Do,  Di  (from  Sect.  2.2),  and  Ci 
are  non-negative  and  that  the  maximum  (absolute)  column  sum  is  always  less  than 


one1.  Equation  (30)  then  results  in  the  forms 

Mo  =  (eo«  ss] 

)  wTPT  (I7_j  -  D0)_1  d0  +  (eo  w)  (I7  -  Ci)-1  uss, 

(31) 

Mi  =  -  (^o  u«s) 

wTPT  (I7_i  -  D0)-2  d0  +  (eo  w)  e $  (I7  -  Ci)-2  Cius*, 

(32) 

M2 

=  (ej uss)  wTPT  (I7_i  +  D0)  (I7_i  -  D0)-3  d0 

+  (©o'w )  eo  (I7  +  Ci)  (I7  —  Ci)  Ciui5, 

(33) 

1With  the  exception  of  degenerate  cases  where  g  (X)  does  not  have  support  over  certain  regions 
of  the  real  line. 
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and 

^3  =  -  (eoUssj  wTPT  (l7_i  +  4D0  +  Dq)  (I7_i  -  D0)_4d0 
+  (e^w)  ^I7  +  4Ci  4-  Of  j  (I7  —  Ci)  4  Cjuss. 

Note  that  Eq.  (31)  should  result  in  juq  =  1. 


(34) 
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Verification  of  analysis 


The  PMF  of  the  starting  time  estimator  derived  in  Sect.  2  is  verified  by  comparing 
the  analytical  distribution  to  one  estimated  through  simulation.  A  Gaussian  shift- 
in-mean  signal  is  considered  where  under  Hq  the  observed  data  are  zero-mean  unit- 
variance  Gaussian  random  variables  and  under  H\  the  mean  changes  to  the  positive 
value  jU.  The  log-likelihood  ratio  is  used  as  the  detector  nonlinearity, 

g(x)=fi(x-£).  (35) 


In  the  simulation  10, 000  trials  are  run  with  a  threshold  h  =  10  where  the  actual 
signal  start  time  is  uniformly  distributed  between  time  samples  101  and  110  to 
ensure  that  the  Page  test  is  in  its  steady  state  before  the  change  occurs.  From 
the  trials,  a  sample  PMF  and  CDF  of  the  error  of  the  starting  time  estimator  of 
Eq.  (2)  are  generated  and  compared  to  that  computed  from  Eq.  (29)  using  100 
quantization  levels  as  shown  respectively  in  Figs.  1  and  2.  Visually  there  is  little 
difference  between  the  estimated  and  theoretical  PMFs  and  CDFs. 

To  further  substantiate  the  validity  of  the  analytical  result  of  Eq.  (29)  and  to  indi¬ 
cate  how  fine  the  quantization  needs  to  be,  the  p-value  of  the  Kolmogorov-Smirnov 
statistic  is  shown  as  a  function  of  the  number  of  quantization  levels  in  Fig.  3.  The  p- 
value  indicates  the  probability  of  exceeding  the  observed  value  of  a  random  variable. 
Thus,  if  the  p- value  is  small,  there  is  a  significant  difference  between  the  observed 
and  estimated  CDF.  As  expected,  the  difference  becomes  significant  as  the  num¬ 
ber  of  quantization  levels  decreases.  For  this  example  (i.e.,  signal  type,  strength 
and  detector  threshold),  the  approximation  becomes  suspect  when  the  number  of 
quantization  steps  falls  below  100.  The  p-value  is  calculated  using  the  asymptotic 
distribution  of  the  Kolmogorov-Smirnov  statistic  as  described  in  Fisz  [8] . 
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Figure  1:  Comparison  of  analytical  (*)  and  estimated  (bar  graph)  PMFs  of  starting 
time  estimator  error. 


Figure  2:  Comparison  of  analytical  (*)  and  estimated  (staircase)  CDFs  of  starting 
time  estimator  error. 
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Figure  3:  p- value  of  the  Kolmogorov-Smirnov  statistic  versus  the  number  of  quanti¬ 
zation  levels  for  Gaussian  shift- in-mean  signal. 


NATO  UNCLASSIFIED 


-  12  - 


SACLANTCEN  SM-295 


NATO  UNCLASSIFIED 


4 

Estimator  performance 


The  performance  of  the  estimator  in  terms  of  the  bias,  standard  deviation  and 
skewness  are  examined  as  a  function  of  signal  strength  for  Gaussian  shift-in-mean 
and  noncentral  chi-squared  signal  types.  The  bias  is  the  first  uncorrected  moment 
(^i),  the  standard  deviation  is, 

c  =  \f^2  -  Ml, 

and  the  skewness  is 

_  /i 3  -  3/X1/Z2  4-  2/if 

H  -  O 

a6 


(36) 

(37) 


4. 1  Matched  detector  nonlinearity 


The  detector  nonlinearity  for  the  Gaussian  shift-in-mean  signal  type  is  the  log- 
likelihood  ratio  as  described  in  Eq.  (35).  For  the  latter  signal  type,  the  data  are 
centrally  chi-squared  distributed  under  Ho  (X*  —  X%)  and  noncentrally  chi-squared 
distributed  under  H\  with  a  noncentrality  parameter  of  6  ( Xi  ~  (6))  and  degrees 
of  freedom  n.  The  locally  optimal  detector  nonlinearity  is  used  with  the  asymptot¬ 
ically  optimal  bias  [9] 

g(x)  =  x  -ts,  (38) 

where 

TJ=n(l  +  2)1°g(l  +  ^).  (39) 

In  both  cases,  the  detector  nonlinearity  is  matched  to  the  actual  signal  strength; 
that  is,  the  actual  signal  strength  is  known  and  used  to  form  the  detector  nonlin¬ 
earity.  Define  the  signal  power  for  the  Gaussian  shift-in-mean  signal  type  as  6  =  fi2 
and  simply  6  for  the  noncentral  chi-squared  signal  type.  The  degrees  of  freedom 
parameter  for  the  noncentral  chi-squared  signal  is  set  to  n  =  2. 

The  bias,  standard  deviation,  and  skewness  of  the  starting  time  estimator  for  both 
signal  types  are  respectively  found  in  Figs.  4,  5,  and  6  as  a  function  of  signal  strength. 
As  expected,  both  the  bias  and  variance  decrease  as  signal  strength  increases.  Ad¬ 
ditionally,  the  bias  for  the  noncentral  chi-squared  signal  is  substantially  larger  than 
the  bias  for  the  Gaussian  shift-in-mean  signal,  particularly  at  lower  signal  strengths. 
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4.2  Mismatched  detector  nonlinearity 


In  most  applications,  the  signal  strength  is  not  known  exactly.  In  these  situations 
a  design  signal  strength  is  chosen  (say  ji  and  6)  and  substituted  into  the  detector 
nonlinearity.  The  bias,  standard  deviation  and  skewness  for  a  design  signal  strength 
fx2  =  S  =  1  are  respectively  found  in  Figs.  4,  5,  and  6  for  the  Gaussian  shift-in-mean 
and  noncentral  chi-squared  signal  types.  Here  it  is  seen  that  the  bias  and  variance 
increase  as  the  signal  strength  decreases  at  a  greater  rate  than  if  the  detector  non¬ 
linearity  were  exactly  matched  to  the  signal  strength,  particularly  for  the  Gaussian 
shift-in-mean  signal. 

When  the  detector  nonlinearity  is  formed  based  on  a  design  signal  strength,  there 
is  a  minimum  detectable  signal  level.  For  the  Page  test,  this  is  usually  chosen  as 
the  signal  strength  that  causes  the  data  transformed  by  the  detector  nonlinearity  to 
have  zero  mean2.  With  a  0  db  design  signal  strength,  this  occurs  at  approximately 
-3  db  for  the  Gaussian  shift-in-mean  signal  and  at  about  -3.6  dB  for  the  noncentral 
chi-squared  signal.  In  these  regions  it  is  clear  that  the  performance  of  the  starting 
time  estimator  is  unacceptable. 

The  bias  and  skewness  both  exhibit  negative  values  for  high  signal  strengths  whereas 
when  the  detector  nonlinearity  is  exactly  matched  to  the  observed  signal  strength 
the  bias  and  skewness  are  always  positive.  The  negative  bias  and  skewness  values 
occur  because  the  probability  of  crossing  the  threshold  immediately  after  the  signal 
starts  without  a  reset  to  zero  is  large.  Thus,  the  starting  time  estimator  will  often 
have  negative  error  values  (which  is  not  necessarily  bad). 


4.3  Detector  threshold 

The  starting  time  estimator  of  Eq.  (2)  is  based  on  the  Page  test  statistic.  Thus,  the 
design  parameters  of  the  Page  test  may  affect  the  performance  of  the  proposed  start¬ 
ing  time  estimator.  The  primary  design  parameter  of  the  Page  test  is  the  threshold 
used  to  declare  signal-presence;  that  is,  the  value  that  the  Page  test  statistic  must 
exceed  for  a  detection  to  occur.  The  detector  threshold  is  usually  chosen  according 
to  constraints  on  the  false  alarm  performance,  with  large  thresholds  corresponding 
to  large  average  times  between  false  alarms.  An  increase  in  the  Page  test  detector 
threshold  results  in  a  decrease  in  the  probability  of  detecting  finite  duration  signals 
and  an  increase  in  the  average  delay  before  detection.  The  effect  of  varying  the 
Page  test  detector  threshold  on  the  bias,  standard  deviation,  and  skewness  of  the 
proposed  starting  time  estimator  for  0  dB  Gaussian  shift-in-mean  and  noncentral 
chi-squared  signals  with  a  matched  detector  nonlinearity  is  found  in  Figs.  7,  8,  and 

2  Because  the  average  number  of  samples  before  detection  becomes  approximately  exponentially 
related  to  the  threshold  for  negative  means  as  opposed  to  linearly  related  for  positive  means. 
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9.  The  only  expected  trend  that  may  be  observed  is  the  increase  in  variance  as  the 
threshold  increases.  Thus,  it  may  be  desirable  to  have  a  lower  threshold  to  reduce 
the  variance  on  the  signal  starting  time  estimator;  however,  the  requirements  on  the 
false  alarm  performance  of  the  detector  may  take  precedence  and  counter  this  desire. 
It  should  also  be  noted  that  for  the  Gaussian  shift-in-mean  signal  the  performance 
does  not  vary  substantially  with  the  threshold. 
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Figure  4:  Bias  for  the  matched  and  mismatched  detector  nonlinearities  for  Gaussian 
shift-in-mean  and  noncentral  chi-squared  signal  types. 


Figure  5:  Standard  deviation  for  the  matched  and  mismatched  detector  nonlineari¬ 
ties  for  Gaussian  shift-in- mean  and  noncentral  chi-squared  signal  types. 
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Signal  Strength  (dB) 

Figure  6:  Skewness  for  the  matched  and  mismatched  detector  nonlinearities  for 
Gaussian  shift-in-mean  and  noncentral  chi-squared  signal  types. 


Detector  Threshold 

Fipre  7:  Bias  for  0  dB  signal  with  the  matched  detector  nonlinearity  for  Gaussian 
shift-in-mean  and  noncentral  chi-squared  signal  types. 
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Detector  Threshold 

Figure  8:  Standard  deviation  for  0  dB  signal  with  the  matched  detector  nonlinearity 
for  Gaussian  shift-in-mean  and  noncentral  chi-squared  signal  types. 


Figure  9:  Skewness  for  0  dB  signal  with  the  matched  detector  nonlinearity  for  Gaus¬ 
sian  shift-in-mean  and  noncentral  chi-squared  signal  types. 
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4.4  Upper  quantile 

The  objective  in  obtaining  an  estimate  of  the  signal  starting  and  stopping  time 
is  usually  to  perform  further  processing  on  the  segment  of  data  that  contains  the 
signal.  It  may  be  desirable  to  ensure,  within  some  probability,  that  the  starting  time 
estimate  is  less  than  or  equal  to  the  actual  starting  time;  that  is,  force  a  negative 
bias  on  the  starting  time  estimator.  For  instance,  the  upper  a-quantile  may  be 
subtracted  from  the  starting  time  estimator  so  that  the  probability  of  observing  a 
positive  error  is  small  (a).  This  provides  the  desired  negative  bias  and  indicates  how 
often  the  starting  time  estimator  misses  the  actual  start  of  the  signal.  The  value  to 
be  subtracted  is  the  positive  integer-valued  quantile  za  where 

a  <  Pr  {M  >  za} 

OO 

=  PM[m] 

m=za 

=  (eo  w)  eo  (I7  —  Ci)_1  Ci°uss.  (40) 

The  value  of  20.01  for  both  signal  types  and  the  matched  and  mismatched  detector 
nonlinearities  is  shown  as  a  function  of  signal  strength  in  Fig.  10.  The  mismatched 
detector  nonlinearities  are  designed  assuming  a  0  dB  signal  strength  and  the  detector 
threshold  was  set  at  h  =  10. 


Figure  10:  20.01  for  Gaussian  shift-in-mean  and  noncentral  chi-squared  signal  types. 


-19- 


NATO  UNCLASSIFIED 


NATO  UNCLASSIFIED 


SaCLANTCEN  SM-295 


5 

Summary 


Many  signal  processing  applications  require  the  detection  of  an  abrupt  change  with 
subsequent  estimation  of  the  actual  time  of  occurrence  of  the  change.  A  starting 
time  estimator  based  on  the  Page  test  statistic  has  been  proposed  for  this  purpose. 
The  PMF  and  first  three  uncorrected  moments  of  the  estimator  error  have  been  de¬ 
termined  analytically  subject  to  a  quantization  of  the  Page  test  statistic  update.  The 
performance  of  the  estimator  for  Gaussian  shift-in-mean  and  noncentral  chi-squared 
signal  types  was  evaluated  revealing  and  quantifying  the  expected  improvement  in 
performance  (i.e.,  decrease  in  bias  and  variance)  as  signal  strength  increased. 

In  many  applications  the  initial  segment  of  the  signal  may  contain  information  cru¬ 
cial  to  the  post-processing.  In  these  situations  it  is  desirable  to  induce  a  negative 
bias  on  the  starting  time  estimator  to  ensure,  within  some  probability,  that  the  esti¬ 
mator  captures  the  start  of  the  signal  (i.e.,  the  estimate  is  less  than  or  equal  to  the 
actual  starting  time  with  some  given  probability) .  Thus,  the  upper  tail  probability 
of  the  error  of  the  starting  time  estimator  was  derived  from  the  PMF  and  used  to 
determine  the  induced  bias  providing  a  given  probability  of  capturing  the  start  of 
the  signal  for  the  Gaussian  shift-in-mean  and  noncentral  chi-squared  signal  types. 
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Abstract 
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